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CONTRACTOR REPORT 


CONFINED SWIRLING JET PREDICTIONS USING A 
MULTIPLE-SCALE TURBULENCE MODEL 

I. INTRODUCTION 


Confined jets are encountered in many types of engineering equipment such as 
coal gasifier, aircraft gas turbine engine, and ramjet combustors. Swirl is often 
introduced to control the mixing of fuel and oxidizer, flame size, and combustion 
intensities. Thus, it is important to be able to predict the flow structure in terms 
of effects of geometric configuration of the combustor, turbulence intensities of inlet 
stream, and swirl strength. 

In the last decade, advances have been made in turbulence modeling such that 
it is possible to predict certain mean and turbulence characteristics of many flows. 
Among the models, the so-called k-e model has been shown to be adequate in a 
variety of flows [1]. However, in complex flows such as confined swirling flows, 
k-e models are known to be deficient [2] due to the neglect of anisotropy in the 
turbulent viscosity and additional turbulence generation arising from streamline curva- 
ture. Several proposals have been raised in the last few years trying to remedy 
this defect. Specifically, a variant in the k-e model is to modify the sink or source 
term in the e-equation by introducing a gradient Richardson number [3] or a flux 
Richardson number [4] to account for the swirl-induced streamline curvature effect. 
Although the detailed interpretation of these two Richardson numbers is different, 
both modifications assume that the turbulence amplification effect of swirl can be 
modeled due to the changes in the large scale of energetic eddies. In other words, 
these modifications assume that there exists a constant energy transfer rate through 
the energy spectrum and that a single velocity and length scale can be used to 
characterize the whole spectrum of turbulence motion. This is thought to be con- 
ceptually weak because the flow field in a confined vortex tube is not in spectral 
equilibrium due to swirl induced destabilizing effects. This also tends to explain 
why these two modifications, and several other modifications developed on similar 
reasoning, only give improvements of the predictions in certain regions of the swirl- 
ing flow [5,6] . 

To circumvent this defect, a so-called "multiple-scale" turbulence model is used 
in the present report to study the confined swirling flows. A basic feature of this 
model is the partition of energy spectrum . This would allow the energy transfer of 
energetic vortexes to be modeled based on the properties of the large scale turbu- 
lence, such as mean velocity gradients and degree of swirl, while the dissipation rate 
is related to the energy transfer by its own action rather than to the overall kinetic 
energy. 

In section II, a brief description of this multiple- scale turbulence model will be 
given and numerical results of a specific confined swirling jet will be presented in 
section III. One major effort of this study is to test the validity of the multiple- 
scale turbulence model by predicting the confined swirling jet. These studies will 
form a basis for multiple-phase closure models. 



II. ANALYSIS 


A. The Mean Governing Equations 

Introducing the gradient-diffusion Boussinesq type approximation into the mean 
Reynolds equation for conservation of momentum, the transport equations can be cast 
into the standard form [7] (cylindrical coordinates): 


a 
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(pU*> + J (rp v +) - yj- 
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( 1 ) 


where <f> is dependent variable, <f> = 1 for mass conservation, and <j> = U,V,W, for 
three velocity components, S^ is the comprehensive source term, and r is the 

phenomenological exchange coefficient. Source terms S^ for these mean quantities 

can be found in Reference 7 and thus are not repeated here. 


B . Turbulence Model 


Turbulence effects which enter the mean governing equations are modeled 
through gradient-type approximation. The eddy viscosity v. is represented by the 

so-called multiple-scale turbulence model [8,9]. This model includes a different set 
of response equations for each of the large-scale energetic eddies and the dissipative 
eddies. The balance equations for the energy and dissipation rate of large-scale 
vortexes and small-scale eddies are derived separately. Specifically, two transport 
equations are provided representing the evolution of turbulent energy associated with 



energy spectrum k^. To close these equations, two transport equations for the 


energy dissipation rate are derived. The dissipation of the energetic eddies, e , can 

P 

be interpreted as the energy transfer rate from the energy-containing eddy motions. 
On the other hand, the dissipation rate of the small eddies, is related to the 


energy transferred by its own action in the high wave number region of the spectrum 
These transport equations are: 
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is the production of energy-containing eddies by the mean gradients and x = k /e . 

^ * 

Note that the time scale used in equation (5) is different from the one used in Refer- 
ence 8. 


This version of multiple- scale model was developed based on the earlier version 
of Hanjalic et al. [8] and has been, tested on several confined recirculating flows [9]. 
The details of the model development can be found in these two papers. The modeled 
coefficients in the turbulence model are established guided by Reference 8 and are 
detailed in Reference 9. These are: 
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Using these model equations, the Reynolds stresses are expressed according to 


u i u j 


h W + 

t[ 2 \3x. 3 x.y 


+ i k 6 ii 


( 6 ) 


where 


k = k. + k 
t P 


and the eddy viscosity is given by 


v. = C (k + k.) 
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with C =0.09 
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(7) 


3 



C. Boundary Conditions 


The experimental work used to evaluate the turbulence model is a confined 
coaxial jet involving sudden expansion with swirl in the annular jet by Roback and 
Johnson [10]. The annular swirling flow was generated by a 30 deg, free vortex 
swirler. The elliptic numerical procedure used for solution requires prescriptions of 
all boundary conditions of the solution domain. Specification of the inlet plane bound- 
ary conditions require special attention because these conditions have a crucial effect 
on the flow field calculated downstream [11-14]. Ideally, inlet conditions should be 
supplied from the experimental data, in this case, at the expansion plane. Unfor- 
tunately the measurements were not performed by Roback and Johnson at this plane 
due to optical difficulties. Calculation starting at the plane upstream of the expansion 
is also unfeasible because the flow field generated from the fluid passing through 
swirl- vanes', is highly complex and three dimensional. Following Leschziner and Rodi 
[11], the computation domain thus starts at the first experimentally measured plane 
which lies just downstream of the expansion at X/D ^ 0.04 (D is the diameter of the 
confined chamber). Distributions of U, V, W and k were available, however e was 
not. Thus, e is estimated from the measured u.’u.' through the standard k-e model 
relationship : 1 •* 


e 
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9 U. 

l 

9 x. 
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( 8 ) 


These two turbulence quantities, k and e, are then used to estimate the inlet profiles 
for multiple-scale turbulence quantities. In this plane the flow is assumed to be 
highly non-equilibrium, thus the non-equilibrium conditions are applied. This implied 
[9] 


and 


k = 


k t 2 k 


= 0.5 


e t 


Along the solid wall, the dependent variables are matched to the usual "wall 
functions" and equivalent local-equilibrium expressions for the turbulence quantities. 
Specifically, the near wall energy transfer rates and are set equal to one 

another. A symmetric axis was specified where the radial velocity and radial gra- 
dients of all other variables were set to zero except near the axis where W was 
deduced from solid body rotation by assuming W to be zero on the axis. 


III. NUMERICAL RESULTS AND DISCUSSION 


Numerical solutions were obtained with a modified version of the TEACH-T 
computer code of Gosman and Iderah [15], The numerical procedure is developed 
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based on finite- volume formulation with hybrid upwind /central-difference scheme for 
convection terms. The computations are made on a 45 x 35 non-uniform mesh with 
concentrations of nodes in the recirculation regions, in the entrance region, and near 
wall regions. Solutions with this mesh were found to be nearly grid-independent [5]. 

When swirl is introduced in the annular flow of a confined coaxial jet, signific- 
ant radial and axial pressure gradients are set up near the central jet exit. At high 
degrees of swirl (swirl number £ 0.5), given the geometry, of the enclosure chamber 
and design of swirl vanes, a central recirculation zone will form. The types of 
recirculation can be of shape characterized by a large central reverse-flow zone (Type 
2) or of shape characterized by an annular zone of recirculation (Type 1) or some 
combination of both (Type 3) depending on the mass flux ratio of annular jet to 
central jet [ 16] . The comparison of the predictions and experiments will be first 
made on the global features of the flow field. 

One of the important parameters to compare in confined flows is the corner 
reattachment length. This length, reported by Roback and Johnson, is approximately 
50 mm which is about one-third the length observed with non-swirling case for the 
same flow condition. Presented in Table 1 is the prediction of this length by multiple- 
scale and single-scale k-e models. Two corrections of the k-e model mentioned in 
Reference 5 are also included. It can be seen that this length is reduced by about 
10 percent between the single-scale and multiple-scale turbulence model and that the 
magnitude of this length is correctly predicted by the multiple- scale turbulence model. 

TABLE 1. COMPARISON OF CORNER RECIRCULATION 

ZONE LENGTH 


Model 

Reattachment Length 
(mm) 

k-e 

v 77 

Launder et al. 

v 80 

Rodi 

v 70 

Multiple-scale 

v 55 

Experiment 

v 50 


Figure 1 shows comparisons of axial velocities along the centerline using k-e 
model and multiple- scale model. It can be seen that the multiple- scale model improves 
the prediction of the strength of the recirculation and predicts a faster recovery of 
the flow than the k-e model downstream of the central recirculation region. Figure 2 
is taken from the previous study [5] of comparison of two k-e model corrections. 
Based on the analysis discussed in Reference 5, the multiple-scale model suppresses 
the diffusion process, thus giving stabilizing effects in the central recirculation 
region. On the other hand, it serves to destabilize the turbulence and speeds up 
the recovery of the centerline velocity to its asymptotic value after the recirculation 
region. The implication is that the multiple-scale turbulence model relates the energy 
transfer rate depending on swirl impartation and flow structure. In the central 
recirculation region, the k-e model tends to connect the energy transfer rate to the 
local mean strain rate too strongly. The multiple- scale model suppresses this ten- 
dency. However, further downstream of the recirculation zone, where the velocity 
gradient is small, since the energy transfer rate of the transfer range eddies does 
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not depend on the mean flow structure [equation (3,5)], multiple-scale model persists 
on the diffusion process, thus enhances the recovery of the flow. 

This behavior can be observed in Figures 3 and 4. Radial profiles of axial 
velocities are plotted at X/D = 0.833 in Figure 3, the multiple- scale model increases 
the size of the central recirculation region and reduces the diffusion. At X/D = 2.5 
(Fig. 4) , augmentation of turbulence by the multiple- scale model smooths the axial 
velocity profile. In fact, the mean flow is almost fully recovered to the fully developed 
pipe flow at this location. Shown in Figure 5 is a comparison of the predicted central 
recirculation zone boundary using the k-e models and the multiple- scale model with 
data of Roback and Johnson. The k-e model, as well as two modifications, produces 
a shorter central recirculation zone. The prediction of the width and the length of 
the recirculation zone is much improved by the multiple-scale model. The recircula- 
tion zone shape, formed in Roback and Johnson's experimental setup, is of type 3 
mentioned above. This type 3 shape is mimicked faithfully by the multiple-scale 
turbulence model. This feature is plotted in Figure 6(a) at X/D = 0.21 and X/D = 

0.42. Also shown in Figure 6(c) are predictions of Sloan [17] of the same case using 
standard k-e, Richardson number correction to k-e, and a non-isotropic algebraic 
stress model, ASM. The multiple-scale obviously has the edge of predicting the 
"three cell" flow structure in the confined swirling flow. 

Figure 7 compares radial profiles of mean tangential velocities at X/D = 0.833 
and 2.5. The single-scale model is seen to have a more rapid decay of swirl compon- 
ent than the multiple-scale model at the upstream location. However, both models 
predict tangential velocity decay to a solid-body rotation prematurely, while the 
experimental data of Roback and Johnson persist in a combined vortex profile even 
at far downstream. The discrepancies between the calculations and measurements 
indicate the deficiency of isotropic eddy viscosity hypothesis used in this level of 
closure. This discrepancy is also shared by the non-isotropic algebraic stress model 
(ASM) prediction [6]. As indicated by Sloan et al., the failure at the downstream 
location probably has to be resolved by solving the full stress transport equations , 
with convective and diffusive terms, rather than the abbreviated ASM. 

Shown in Figure 8 are predictions of turbulence intensity at several downstream 
locations. Both models follow the trend of experimental data fairly well except that 
the peak values are underpredicted. Multiple-scale model seems to have a more 
realistic presentation of the fluctuation near the wall region. Again, the higher level 
of turbulence intensities at the downstream region is that the kinetic energy contained 
by the transfer eddies (or small-scale eddies), k t> is independent of the mean strain 

rate. Although the production of k is very small here, the production of k. still 
persists. P 

Recently the sensitivity of the numerical prediction to the inlet boundary condi- 
tions has gained increased concern [11-14]. It is instructional to examine the 
response of the flow field to variations in turbulence quantities of inlet streams. 

Since the mean variable profiles are taken from experimental data, the uncertainty 
associated with mean field inlet boundary conditions would be minimal. The specifica- 
tion of turbulence quantities, such as k^, k^, Cp and e^., on the other hand, is in 

doubt. Computations compared above were performed using "non-eqilibrium" spectral 
shape. In order to examine the effect of the assumed spectral shape of turbulence, 
calculations were made with "partially-eqilibrium" and "equilibrium" inlet spectra. In 

the partially-equilibrium spectrum, the kinetic energy is separated such that k = 

P 

0.7 k and k { = 0.3 k. In the equilibrium case, both inlet Cp and e t were set equal 
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to e and the division of the energy spectrum is moved to high wave numbers such 
that only a very small portion of the turbulence kinetic energy is contained in k^.. 

As seen from Figures 9 and 10, the difference of mean flow field predictions between 
using non-equilibrium and partially-equilibrium initial spectra is not appreciable. The 
prediction starting with equilibrium spectral shape, however, shows substantial depar- 
ture from that of non-equilibrium shape. The kinetic energy spectrum in most inlet 
planes of complex flows cannot be equilibrium. Thus, it is recommended that non- 
equilibrium or close to non- equilibrium spectral shape should be used to start the 
calculation . 


IV. CONCLUDING REMARKS 


Despite the fact that the proficiency of turbulence model predictions may be 
obscured by the influences of boundary conditions , numerical diffusion , grid system 
and oscillatory phenomena, significant benefits can be gained in elucidating the rela- 
tive advantages of various turbulence model predictions. It is demonstrated in this 
study that the multiple-scale turbulence model exhibits great potential in swirling 
flow predictive capability. The preliminary success of the multiple-scale turbulence 
model stems from recognizing that the turbulence in swirling flow is highly non- 
spectral-equilibrium, and, that the different energy transfer rates for different 
eddies should be modeled separately. However, the multiple- scale turbulence model 
still adopts the isotropic eddy viscosity formulation. To take into account the 
anisotropy of the flow turbulence, full Reynolds Stress Model has to be used. Incor- 
poration of the multiple- scale turbulence model, which takes into account the non- 
equilibrium spectral energy transfer rate and simplified Reynolds Stress Model, should 
pave a promising avenue for numerical modeling of swirling flows. 

One last note is that although much research has to be aimed at improving 
predictive capabilities of current turbulence models, the development of the numerical 
efficiencies and accurate differencing schemes for validating and optimizing turbulence 
models should not be overlooked. More recent works [18-20] have been devoted to 
these developments. Also, fundamental experimental investigations including the 
detailed inlet boundary conditions measurements and characterization of the discharg- 
ing flow from swirlers remain an absolute necessity for guiding and testing computa- 
tional models. 
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